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1 .  Purpose  of  Study 

The  main  objective  of  this  study,  as  originally  set  forth  in  work 
statement  79-277,  is  to  define  quantitatively  the  range  and  accuracy  of 
mu Itlspectra I  lidar  techniques  for  remotely  sensing  the  visibility,  sire 
distribution,  concentration,  and  liquid  water  content  of  various  types  of 
fogs.  As  outlined  in  the  above  work  statement,  the  study  was  to  Include 
three  Ildar  wavelengths  (0.53  pm,  1.06  pm,  and  10.6  pm).  Radiaticr  and 
advectlon  fogs  with  visibilities  of  30  —  1 000 meters  were  to  be  investi¬ 
gated  to  ascertain,  first,  from  a  single  wavelength,  how  well  the  above 
parameters  could  be  determined,  and  second,  using  information  frcm  three 
wavelengths,  whether  the  errors  in  the  determinations  could  be  reduced. 
Finally,  the  maximum  useful  ranges  for  1.06  pm  and  10.6  pm  wavelengths 
were  to  be  found  for  the  given  range  of  visibilities,  and  the  importance 
of  multiple  scattering  was  to  be  determined. 

The  original  work  order  was  subsequently  modified  (personal  communi¬ 
cation  with  Dr.  J.  S.  Randhawa),  and  the  study  was  narrowed  down  to  the 
following  objectives: 

a.  Four  lidar  wavelengths,  representing  currently  available  systems 
within  the  wavelength  region  of  i nterest,  were  to  be  studied.  The  wave¬ 
lengths  are  0..53  pm,  1.06  pm,  3.8  pm,  and  10.6  pm. 

b.  Lidar  performances  and  s igna I -to-noi se  (S/N)  ratios  are  to  be 
ascertained  for  representative  systems  at  these  wavelengths. 

c.  Maximum  useful  ranges  at  the  4  wavelengths  were  to  be  determined. 

d.  Methods  to  obtain  extinction,  back-scatter  and  total  liquid  water 
content  were  to  be  investigated. 

It  was  agreed  that  the  question  of  obtaining  size  distributions  and  the 
effects  of  multiple  scattering  would  not  be  considered  here,  due  to  the 
complexity  of  these  problems  and  the  time  restrictions  on  the  present  work. 


2.  Discussion  of  Fog  Models 

In  order  to  perform  the  required  calculations  for  this  study,  it  is 
necessary  to  model  the  distribution  of  droplet  sizes  composing  the  fogs, 
as  the  optical  parameters  (t.e.,  extinction,  backscatter,  optical  depth) 
wi 1  I  depend  upon  the  chosen  distributions.  While  many  fog  distributions 
are  described  in  the  literature,  it  was  decided  that  those  given  by 
Pinnick  et  at.  (!978)  were  most  relevant  and  convenient  fcr  +ho  r „r poses 
of  this  study.  pi nnick  et  at.  have  conveniently  fit  observe i  r i zo  distri¬ 
butions  to  a  superposition  of  two  log-normal  distribution  functions  i  '  c 
Pig.  I)  and  have  listed  the  values  of  the  apprcpr iate  parameters  fcr  ■  -'h. 
The  advantage  of  the  Pinnick  distributions  is  that  the  analytic  form; 
renders  them  easily  amenable  to  computer  calculations.  rurthc-rmcre,  by 
varying  the  parameters  of  the  log-normal  functions,  variations  of  tl,e 
initial  distributions  may  readily  be  generated  with  ccrr '  de-abl  y  '•"■•■■s 
computer  effort. 

The  log-nornal  distribution  function  is  given  by 

d?T— -  =  - L- —  exp  C -(.In  r  -  In  r)2/2a2]  (  I  ) 

dr  (2»)V 

Here,  n(r)  is  the  number  of  particles  per  unit  volume  of  radius  r,  r  is  the 
mean  radius,  and  a  is  the  natural  log  of  the  standard  deviation  of  r. 

The  first  distribution  given  by  Pinnick  et  at.,  which  we  here  refer  1o 
as  distribution  A  or  dist.  A,  Is  a  superposition  of  two  log-normal  functions 
with  parameters  r  =  0.3  pm,  o  =  1.87,  and  total  liquid  water  content  (LWC) 

=  I.65xl03  pg/m3,  and  r  =  4.4  pm,  o  =  1.47,  and  LWC  =  I.60xl05  pg/m  .  i he 
second  distribution,  distribution  B  or  dist.  B,  has  parameters  r  =  0.6  pm, 
o  =  2.5,  and  LWC  =  7.52x|04  pg/m3  and  r  =  5.0  pm,  o  =  1.65,  and  LWC  = 
4.37x|03  yg/m3.  Plots  of  distributions  A  and  B  are  shown  in  Fig.  I. 
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Figure  1.  Size  distributions  (from  Pinnick  et  al  1978)  used  for  the 
calculations  presented  in  this  report.  (See  text  for  a 
discussion  of  the  modifications  to  the  above  two  distribu¬ 
tions  that  were  also  employed.) 


Besides  the  standard  distributions  A  and  B  as  listed  above,  new  distri¬ 
butions  were  generated  as  follows.  The  liquid  water  content  of  the  two 
individual  log-normal  functions  raking  up  distribution  A  are  In  the  rat io 
of  about  97:1,  while  those  making  up  distribution  B  are  in  the  ratio  5.6: I . 
By  keeping  the  same  the  Individual  normalized  log-normal  functions  which 
make  up  each  distribution,  but  varying  the  ratio  of  the  LWC's,  a  wide  range 
of  new  size-distribution  functions  may  be  generated.  ~h:s  method  f  os  I  he 
advantage  of  not  requiring  repeated  calculations  for  each  row  d i si r i but  ion . 
Once  the  optical  properties,  such  as  volume  extinction  or  b acksc •'c 
coefficient,  have  been  determined  for  each  individual  log-normal  function 
in  each  distribution,  these  properties  need  be  simply  multiplied  by  the 
proper  constants  to  provide  tne  required  ratios  of  the  LWC's.  Liquid  water 
content  ratios  of  the  two  log-normal  functions  making  up  rjCh  oistributlo'’ 
were  varied  from.  1000:1  to  0.1:1.  A  few  of  the  resulting  size  d i s 1 r i ration 
functions  are  shown  in  Figs.  2  through  7. 

3.  Computation  of  Optical  Properties 

Use  of  the  I idar  equation  requires  knowledge  of  the  volume  extinction 

coefficient,  Q^y>  and  the  volume  backscatter  coefficient,  g,  in  order  to 

perform  calculations  of  the  returned  signal  and  the  S/N  ratios  to  be 

expected.  The  calculations  of  and  g  were  performed  as  follows. 

Calculations  of  the  backscatter  cross  section,  g(r),  and  extinction 

cross  section,  Q^(r) ,  for  a  particle  of  radius  r  were  made,  using  the 

standard  Mie  equations  and  the  appropriate  complex  indices  of  refracticn. 

These  are  related  to  the  volume  coefficients  through  the  relationships 

00 

e  =  f  c  (r)  d-”-(£l  dr 
hTV  J  wT  dr 


(2a 
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Figure  4.  Log-normal  distribution  function  as  calculated  for  various 
sizes  (for  Distribution  A,  mass  ratio  =  100/1). 
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Figure  7.  Same  as  in  figure  6  (Distribution  B). 
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and 


f 

f 

f 


6(r>  dr  , 

dr  ' 


where  [d«(r)H/dr  is  one  of  the  Pinnlck  distributions  or  Its  modifications, 
as  discussed  in  section  2.  In  practice,  the  integrations  were  carried  out 
numerically  for  each  of  the  two  log-normal  distributions  in  dist.  A  and 
dist.  B.  Then  the  total  LWC  for  each  log-normal  distribution  was  conouted 


from  the  equation 

LWC 


, ,,  3  dn(r) 

=  p  4/3  irr  — - - 

dr 


dr 


where  p  is  the  density  of  liquid  water  (=1.0).  Then  the  individual  log¬ 
normal  distributions  for  the  two  fog  types  (dist.  A  and  dist.  B)  were  multi¬ 
plied  by  the  proper  constants  to  bring  the  individual  LWC’s  into  the 
required  ratios.  These  constants  were  then  used  to  combine  the  values  of 
Q^v  and  8  for  the  individual  log-normals,  to  give  the  required  values  for 
dist.  A  and  dist.  B.  Once  the  log-normals  were  combined,  they  were  then 
renormalized  so  that  dist.  A  and  dist.  B  each  represented  a  total  liquid 
water  content  of  I  gm/m3  for  the  modified  Pinnick  distributions.  The 
liquid  water  content  for  the  original  Pinnick  distributions  was  set  as 
given  by  Pinnick  et  at.,  i.e.,  0.161  gm/m3  for  dist.  A  and  0.512  gm/m3  for 
dist.  B.  Thus  the  modified  Pinnick  distributions  represent  approximately 
2-5  times  the  LWC  as  originally  measured  by  Pinnick.  However,  these 
modified  distributions  can  be  used  for  any  desired  LWC  by  simply  multiply¬ 
ing  the  present  results  by  the  proper  renormalizing  constant.  The  calcu¬ 
lations  of  S/N  ratios  presented  in  section  3  used  the  original  LWC  as 
given  by  Pinnick  et  at. 

Figures  8  and  9  show  for  a  fixed  LWC  =  I  gm/m3,  as  a  function  of 

the  various  modified  Pinnick  distributions  A  and  B  for  each  of  the  four 
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MODIFIED  DISTRIBUTION  A 


selected  wavelengths.  The  abscissa  in  these  figures  represents  He  "  :ss 
ratios  cf  the  two  log-normal  functions  making  up  each  d i st r i  ‘  ut icr . 

These  figures  are  Interesting  in  that  they  show  clearly  the  ir: or tance 
of  the  drop-size  spectrum  In  determining  the  resulting  optical  preper+los 
cf  the  fog,  particularly  for  the  modifications  of  distribution  A.  re,  .it 
the  two  shorter  wavelengths  (X  -  0.55  pm  and  X  =  I .00  pm),  the  extirctior 
changes  by  about  one  order  of  naanl+ude  over  ■‘•he  range  of  modi  Med  distri¬ 
butions  considered.  It  is  important  tc  ran  erber  tha^  all  diet  r  ib-jl  is  ns  in 
these  figures  have  been  normalized  to  the  same  I.WC,  so  the  yf‘oc's  are  \.r 
solely  to  the  different  drop-size  spectra.  For  d i str i but  ion  3,  the  c’-.r  oos 
are  loss  dramatic,  being  approximately  a  factor  cf  2  over  the  complete  set 
of  size  distributions  at  the  three  shorter  wavelengths.  Of  perhaps  core  sig¬ 
nificance  is  the  fact  that,  at  10.6  pm,  and  to  so" o  extent  at  '■  .0  p",  for  all 
modifications  of  both  dist.  A  and  dist.  B,  the  volume  extinction  is  nearly 
constant.  This  suggests  the  possibility  of  using  the  determination  of  Q,.,r  as 
a  direct  measure  of  the  LWC  o*  the  fog.  This  :s  discussed  in  section  5. 


4.  Lidar  Techniques 

a.  Lidar  Equation 

The  lidar  equation  for  a  pulsed  coaxial  system  can  be  written  as 

c<7. 

P  ~  A  — 5-  B  exp  -  2 
r  2B2 


(3) 


where 


A 

R 

6 


and 


TV 


o 

the  radiation  backsca ftered  to  the  lidar  receiver, 

the  area  of  the  receiver, 

the  energy  In  the  transmitted  laser  pulse, 

the  range  to  the  scattering  pulse, 

the  differential  backscatter i ng  cross  section  per  unit 
volume  of  scattering  medium, 

the  medium  extinction  cross  section  per  unit  scattering  volume. 
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Three  assumptions  have  been  made  In  the  development  of  this  equation. 
The  first  of  these  is  that  the  scattering  process  Is  the  result  of  single 
scattering;  multiple  scattering  Is  Ignored.  The  second  assumption  Is 
that  molecular  scattering  can  be  Ignored  when  ccrpared  to  the  fog  droplet 
scatter.  This  follows,  since  the  differential  molecular  volume  backscattcr 
cross  section  at  sea  level  Is  approximately  10  m  /sr  compared  to 
10  3  m  Vsr  for  fog.  These  values  are  for  a  wavelength  of  0.55  p,  wk!ch 
would  magni'y  the  molecular  scattering  contribution  when  compared  with 
values  at  longer  wavelengths. 

The  third  assumption  Is  that  the  attenuation  of  the  laser  pulse  is 
brought  about  solely  through  fog  droplet  scattering  and  absorption  losses. 
The  laser  wavelengths  have  been  chosen  such  that  molecular  absorption  can 
be  ignored.  Molecular  scattering  can  be  neglected  using  arguments  similar 
to  those  used  in  discussing  the  volume  scattering  cross  sections. 


b.  Lidar  Parameters 

Four  lidar  systems  have  been  chosen  for  analysis;  these  lidars  c cerate 
in  spectral  regions  of  interest  to  the  sponsoring  agency.  The  parameters 
listed  below  for  the  lidars  represent  available  components  near  the  state 
of  the  art.  Wavelengths  are  listed  in  microns,  and  energy  is  measured  in 
jou I es . 

TABLE  I 


A(p) 

Energy 

U) 

System 

Transmi ttance 

FI  Iter  Width 

<y> 

Detector 

.53 

.2 

.6 

io"3 

S-20 

1.06 

1 

.4 

io"3 

Avalanche  Diode 

3.8 

1 

.6 

2* 1 0-2 

1  nSb 

10.6 

10 

.8 

2x|0-2 

HgCdTe 

i 


j 


In  each  case,  the  receiver  aperture  was  taken  to  be  0.5  m  In  diameter, 
the  receiver  field  of  view  to  be  10  *  sr,  and  the  signal  processing  tand- 
wldth  to  be  I07  Hz. 

c.  Lidar  Performance 

Signal-to-nolse  (S/N)  calculations  have  been  performed  for  the  lidar 
system  described  above.  It  was  assumed  In  these  calculations  +kat 
Ildars  were  immersed  in  a  ground  fog  of  the  type  described  in  Fie.  1, 
taken  from  the  paper  by  Pinnick  et  al .  ('978).  In  order  to  ■.'."side-  a 
wider  range  of  fog  conditions,  the  analysis  was  also  perfor'-ed  'or  #ive 
fogs  having  the  same  size  distribution  but  with  liquid  water  consents 
ranging  from  a  factor  of  two  larger  than  the  original  distribution  Jo 
a  factor  of  1/8  the  original  fog. 

The  differential  volume  backseat ter tng  cress  sections,  8  (m-1sr-1), 
and  the  extinction  coefficients,  Q ^  (m-1)  for  the  basic  distributions 
are  as  fol lows. 

TABLE  II 

(See  Figure  I  for  description  of  distributions  A  and  B) 


Qmy 

Wavelength  - 

(Dist.  A)  (Dist.  B) 


6 


(Dist.  A)  (Dist.  B) 


.53  .443E-0I 

t.06  .449E-0I 

3.8  .5I9E-0I 

10.60  .237E-0I 


.  l04Et00 
.  !08EtOO 
. I I8E+00 
.695E-0I 


. 206E-02 
. 200E-02 
.S78E-03 
.2I9E-04 


. 526E-02 
. 499E-02 
. 244E-02 
. 497E-04 


In  order  to  determine  the  S/N  ratio,  it  was  necessary  to  make  an  assump¬ 
tion  as  to  the  levels  of  the  background  radiation.  These  values  were  deter¬ 
mined,  in  the  case  of  those  Ildars  working  at  wavelengths  4  microns  and  less. 
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by  scattering  the  appropriate  solar  Irradlance  isotropically.  These  values 
agree  within  a  factor  of  2  of  those  suggested  in  Wolfe  and  Zissis  (1976). 

It  was  assumed  for  the  10.6  micron  system  that  the  background  level 
originated  from  the  fog  acting  as  a  black  body. 

The  signal  level  as  a  function  of  range  was  computed  using  Eq.  (3)  and 
the  system  parameters  listed  in  Tables  I  and  II  for  each  of  the  two  basic 
size  d  i  str i but  ions .  The  system  noise  was  determined,  according  to  the  type 
detector  used  in  each  system,  from  typical  manufacturer  spec i f icat ions . 

The  system  noise  is  somewhat  range-dependent,  since  it  depends  upon  both 
the  constant  intrinsic  detector  noise  and  background  noise,  as  well  as  the 
time-varying  signal  shot  noise.  The  detection  process  was  considered  +o 
be  incoherent,  and  all  the  results  are  valid  for  single  pulses  (t.e.,  no 
S/N  enhancement  through  averaging). 

The  results  of  these  calculations  are  presented  in  Figs.  10  through  17. 

5.  Data  Evaluation 

a.  Determination  of  Optical  Properties  of  Fog 

Reference  to  the  curves  of  s igna I -to-no i se  ratio  as  a  function  of  fog 
penetration  indicates  that  useful  Ildar  signals  can  be  measured  only  over 
distances  of  the  order  of  a  few  hundred  meters.  It  should  also  be  recog¬ 
nized  that  the  fog  droplet  concentration  as  well  as  size  distribution 
will  fluctuate,  both  with  time  and  with  position  within  the  fog.  Con¬ 
sequently,  It  seems  reasonable  that  effective  values  of  both  6  and  Q^y 
be  derived  from  the  Ildar  measurements  which  represent  both  temporal  and 
spatial  averages  of  these  cross  sections. 

Equation  (3)  can  be  written  explicitly  In  terms  of  the  spatial  and 
temporal  fluctuations  in  8  and  Q^y  as  follows. 
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The  time  average  P^{F^,t)  will  reduce  the  term  con+air.ina  Q'{F  ,t)  tc;  n 
nenllaiblr-  value  under  the  assumption  that  the  spatial  average  of  the  random 
variable  Q'(R,t )  performed  In  the  exponential  integral  Is  suitably  small. 
cnuation  (4)  then  takes  the  form 
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Writing  it  logarithmically  yields 


a 


where  C  =  ACJ./2 ,  and  the  averages  are  time  averages  of  the  returns,  all  from 
range  A  .  We  will  drop  the  bars  in  the  subsequent  discussion.  Now,  if  the 
value  of  can  be  considered  to  be  a  mean  over  the  distance  from  0  to 

,  it  can  be  taken  out  of  the  integral,  making  the  last  term  2  .<y,  {Fy)  , 
where  is  the  effective  average  volume  extinction  over  the  path  from 

0  to  F  . 

n 

Now,  we  assume  we  have  returns  from  a  large  number  of  distances  F.  ,  n  = 

n 

1,  2,  ...,  /V.  We  can  then  compute  the  effective  averages  for  3  and  .4’  from 
the  entire  set  of  returns  in  a  least-squares  sense  in  the  following  manner. 

We  take  the  sum  of  all  the  returns  for  all  R  ,  n  =  1,  2,  N.  Thus 

n 

l  \ln  PR(Rn)  =  In  C  -  Z  In  Rn+  In  {$(*„)  -  2  Q^R  )  *1 . 
n=\  —I 
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We  now  ask  the  question  as  to  what  is  the  best  single  value  for  B  and  y.,.,,  to 
fit  all  the  measurements  in  a  least-squares  sense.  Thus, 
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n  =  i 


l  n  C  +  2  l  n  R 
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In  i'  +  2 


where  we  want  to  minimize  A.  Therefore,  take  3A/36  and  3A/:)..V.  and  set  each 
to  zero,  then  solve  the  pair  of  equations  for  B  and  to  yield  the  best  fit. 
This  is  straightforward  and  yields 
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These  two  expressions  should  yield  values  of  In  6  and  which  fit  the  mea¬ 
surements  of  Pn  at  various  ranges  Ft  ,  best  in  a  least-squares  sense.  Note 

R  n 

that  they  fit  the  logarithmic  form  of  the  lidar  equation  best. 

Now,  one  could  get  "best-fit"  values  of  $  and  Q  over  different  segments 

TV 

of  the  entire  path.  Thus,  the  entire  path  could  be  broken  up  into  segments 
RN ),  ...,  etc.  Over  each  segment,  say  a  best  fit  for  Q^y  and 
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l>i£  could  be  obtained.  The  value  for  CfV  could  then  be  used  to  get  the 
incident  beam  strength  on  the  next  segment,  PN2 ,  and  a  similar  procedure  could 
be  used  to  compute  best-fit  values  over  the  next  segment.  This  could  be 
repeated  as  required. 

b.  Detcvrrinaticn  of  Liquid  Water  Content 

r  toure  18  demonstrates  more  clearly  the  relationship  o.  -  <cuid  .■>  r 
content  to  0TV  of  the  fog.  In  this  figure,  the  values  of  <?...••,  for  o' I  of 

the  mod i f i cat i ons  of  distributions  A  and  B  are  plotted  fer  ouch  of  four 

wavelengths.  Noting  that  all  of  the  calculations  of  Q  were  for  a  fixed 
LWC  of  I  gm/m3,  it  is  clear  that,  in  order  for  a  determination  of  Q  to  be 
Interpreted  as  a  measure  of  +he  LWC  of  the  fog,  there  must  be  a  unique,  or 
nearly  unique,  relationship  of  Q ^  to  LWC.  This  relationship  must  be  inde¬ 
pendent  of  the  unknown  size  distribution.  This  figure  shows  that,  at  1 ne 
two  longer  wavelengths,  and  particularly  at  the  longest  wavelength  (10.6  ym) , 
the  relationship  is  indeed  nearly  unique.  Thus,  for  the  present  set  of  cal¬ 
culations,  all  size  d  i  str  i  but  ions  employed  yielded  a  value  of  -  -18  m 

for  an  LWC  =  1.0  g/m3.  If  the  LWC  were  doubled,  the  resulting  QfV  would  also 
be  doubled.  It  is,  of  course,  possible  that  real  fogs  may,  in  fact,  have 
size  distributions  so  grossly  different  from  the  wide  range  of  possibilities 
considered  here  that  this  apparent  relationship  would  break  down.  However, 
from  the  broad  range  of  possibilities  considered  here,  it  does  appear  that 
a  measure  of  10.6  ym  will  provide  a  reasonable  estimate  of  the  LWC 

in  the  fog.  These  results  are  consistent  with  those  of  Pinnick  et  at. 

(1978b)  and  with  the  theoretical  study  of  Chylek  (1978). 
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6.  Conclusions 


The  result's  of  This  study  can  be  summarized  by  the  following  ■, !  at  >  rents, : 

a.  The  penetration  of  the  I i dar  pulse  into  the  fog  is  limited  to  a 
few  hundred  meters  due  to  the  severe  attenuation  of  the  pulse  by  flu.  foe. 

The  penetration  depth  increases  with  increasing  wavelength,  as  shown  in 
Figures  10  through  17,  in  which  the  usefulness  of  the  returned  Mgnal  : s 
fasured  by  its  s  igna  I -to-nci  se  ratio. 

b.  Methods  are  suggested  to  obtain  the  averages  of  £  -v'c  w".ich 
utilize  both  multiple-pulse  temporal  and  least-squares  spa!  id  i  averag  i  r.p 

p  rocedures. 

c.  It  has  been  shown,  consistent  with  the  literature,  that  llnuio 

water  content  can  be  effectively  measured  from  a  knowledge  of  u-  ing 

u  ► 

a  10.6  p  Ildar.  Further,  these  studies  indicate  that  I  i dar  operating  at 
a  wavelength  of  4  p  would  also  be  useful  for  this  purpose. 

7.  Recommendations  for  Future  Work 

a.  It  is  recommended  that  experimental  verification  of  the  techniques 
for  obtaining  Qj,y,  £,  and  IWC,  as  outlined  in  this  report,  be  underfaken  to 
determine  the  practical  feasibility  of  our  results.  The  required  laser 
systems  and  components  assumed  in  this  study  are  currently  available,  and 
therefore  such  a  program  Is  now  possible. 

b.  At  optical  depths  exceeding  a  few  tenths,  which  in  most  fons  is 
reached  within  tens  of  meters  or  less,  multiple  scatter  effects  may  become 
important.  It  is  recommended  that  theoretical  and  experimental  studies  of 
these  effects  be  undertaken.  The  theoretical  treatment  of  multiple  scatter 
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of  laser  beams  is  quite  complex,  although  some  models  have  been  de\ eloped. 
Nevertheless,  this  problem  cannot  be  adequately  studied  in  a  short  time  period 
such  as  was  available  for  the  current  work,  and  it  is  recommended  that  adt  ;uate 
financing  and  time  be  made  available  for  a  proper  treatment  of  this  problem. 

c.  The  values  of  Q ^  and  8  at  the  four  wavelengths  used  in  this  study 
contain  some  information  as  to  the  size  distribution  of  the  droplets  composing 
the  fog.  A  complete  set  of  these  parameters  would  represent  8  pieces  of  infor¬ 
mation,  in  the  absence  of  errors  of  any  type.  The  presence  of  errors  (c.g., 
measurement  errors,  numerical  quadrature  errors,  etc.)  will  undoubtedly  reduce 
the  amount  of  independent  data  below  the  original  8  pieces.  Nevertheless,  it 
would  appear  that  at  least  an  estimate  of  the  size  distribution  could  be 
obtained  through  appropriate  inversion  techniques.  It  is  recommended  that 
this  problem  be  studied  to  gain  an  insight  as  to  how  well  the  size  distribution 
functions  may  be  determined. 
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